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We present a detailed investigation of the acene series using high-level wavefunction theory. Our 
ab-initio Density Matrix Renormahzation Group algorithm has enabled us to carry out Complete 
Active Space calculations on the acenes from napthalene to dodecacene correlating the full tt- valence 
space. While we find that the ground-state is a singlet for all chain- lengths, examination of several 
measures of radical character, including the natural orbitals, effective number of unpaired electrons, 
and various correlation functions, suggests that the longer acene ground-states are polyradical in 
nature. 
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I. INTRODUCTION 

The acenes are the series of ladder-like compounds con- 
sisting of linearly- fused benzene rings (Fig. [I]) [U, ■ Due 
to their technological potential [1, U@ and their intrinsic 
value as models for more complex conjugated molecules, 
they have been the subject of many theoretical and ex- 
perimental investigations [1, 0, S, B dill!, [H [H, [3 ■ 
In a number of recent studies, it has been proposed that 
longer acenes may possess an unusual electronic ground- 
state that is not the simple closed-shell singlet suggested 
by molecular orbital arguments. Based on extrapolating 
the experimental singlet-triplet gap of the acenes up to 
pentacene, Angliker et al. Q predicted that the ground 
state of higher acenes from nonacene upwards would be 
a triplet. Density functional calculations by Houk et al. 
[ni also predicted a singlet-triplet cross-over. However, 
Bendikov et al. [T5| noted that the restricted singlet den- 
sity functional ground-state would become unstable to 
an open-shell singlet, or singlet diradical, configuration 
for acenes longer than hexacene. The open-shell singlet- 
triplet gap for the longest acene studied (decacene) was 
estimated as ranging from 1.5 (BLYP/6-31G(d)) to 5.7 
(B3LYP/6-31G(d)) kcal/mol. 

Despite these intriguing findings, the density func- 
tional results leave many interesting questions unan- 
swered. For example, how diradicaloid are the acenes 
really compared to conventional diradical systems? As 
we go to longer acenes, might we expect to find tri- and 
even higher polyradical ground states? And if so, how 
do we understand the electronic structure and bonding 




FIG. 1: The first few members and the unit cell of the acene 



in these states? Such questions, which probe the essential 
many- electron character of di- (and indeed poly-) radical- 
ism, are not easily answered through density functional 
theory based on a single Kohn-Sham determinant. 

For this reason we have decided to explore the nature 
of the acene ground-state using high-level wavefunction- 
based electronic structure theory. The many-electron 
correlations in radical wavefunctions tell us about the 
coupled simultaneous motions of the electrons. Con- 
ceptually, singlet states with unpaired electrons require 
multi-configurational wavefunctions [Tsl [H, [13, [H, [H, 
[20l . [21I [2^. [23| as used, for example, in the CAS (com- 
plete active space) family of methods [24|. In the acenes, 
the ideal choice of active space would be the complete 
TT-valence space, i.e. the set of all conjugated pz orbitals. 
However, the exponential cost of traditional CAS meth- 
ods as a function of the number of correlated orbitals 
and electrons renders calculations with the complete tt- 
valence space impossible for acenes much longer than 
napthalene, which already has 10 conjugated orbitals and 
electrons. Consequently, earlier CAS calculations could 
only use an incomplete tt- valence space [H, . 

The Den sity Matrix Renormahzation Group (DMRG) 
algorithm ^Q, [23, [ll] provides a way to overcome the 
traditional exponential complexity of CAS methods in 
long molecules such as the acenes We have re- 

cently developed a local ab-initio DMRG method that 
computes an essentially exact CAS wavefunction with an 
effort that scales only quadratically with the length of 
the molecule (3T| . Consequently, we can now extend the 
range of traditional CAS calculations in long molecules 
from about 10 orbitals to 100 active orbitals or more. 
In the current work we apply our ab-initio DMRG al- 
gorithm to the acene series from napthalene (2-acene) to 
dodecacene (1 2-acene), in all cases correlating exactly the 
complete 7r-valence space. First we revisit the question 
of the relative stabilities of the singlet and triplet states. 
Then, using our correlated wavefunctions, we embark on 
a detailed study of the radical nature of these systems. 
We find intriguingly that the higher acenes are not only 
diradicals, but possess increasing polyradical character. 
By explicit visualisation of the electron correlation, we 
uncover a coupled motion of the electrons that gives a 
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new picture of bonding in molecules with extended con- 
jugation, showing that even systems such as the acenes 
can continue to provide fertile sources of surprising elec- 
tronic structure. 



II. COMPUTATIONAL METHODOLOGY 

In the present context, we can regard the DMRG as an 
efficient way to exactly correlate, in the sense of full con- 
figuration interaction, the electrons in the active space. 
Details of the DMRG algorithm as implemented in our 
Block code are given in Ref. [3l|, [H, [11]. Active 
space full configuration interaction theory is sometimes 
referred to as CASCI. Recall that CASCI is the same as 
the more common CASSCF (complete active space self- 
consistent field) method [131 but lacks the step of orbital 
optimisation. Orbital optimisation is possible within the 
DMRG but has not been used here. For the molecules 
in this work, the DMRG energies are converged to better 
than 0.1 kcal/mol and would be identical to the so-called 
CASCI energies if it were possible to compute these in 
the traditional manner. 

DMRG calculations on the acenes were performed at 
the UB3LYP/6-31G(d) P, [H optimised singlet and 
triplet geometries which were essentially the same as 
those used by Bendikov et al. [l2| (see supporting in- 
formation for details). These structures have D2/1 point- 
group symmetry. The rung C-C bonds are somewhat 
longer than the ladder C-C bonds and the ladder C-C 
bonds display increasing bond alternation towards the 
ends of the chain. For example in singlet decacene, the 
rung and ladder C-C bonds were 1.464A and 1.405A re- 
spectively at the middle of the chain, while the difference 
in successive ladder C-C bonds lengths was 0.058A at the 
end of the chain as compared to O.OlOA at the middle. 

The active space was chosen to be the complete tt va- 
lence space, consisting of all conjugated carbon pz or- 
bitals, and all tt electrons were correlated. The a elec- 
trons were treated within a frozen-core approximation 
using the restricted Hartree-Fock orbitals. The calcula- 
tions used either the minimal ST0-3G basis 30| (up to 
dodecacene) or Dunning's double-^ DZ [13, [38| basis (up 
to hexacene) as indicated. In the case of the DZ basis, 
two Pz orbitals per carbon were used to make a "dou- 
ble" complete 7r-valence space. Thus whereas e.g. the 
DMRG/ST0-3G calculations for pentacene correspond 
to a (22, 22) CASCI, the DMRG/DZ calculations would 
correspond to a (22, 44) CASCI. ^811] 



31G(d) and UBLYP/6-31G(d) singlet-triplet gaps (using 
open-shell wavefunctions where stable) as first reported 
by Bendikov et al. [T^ . which we have recomputed and 
extended to the complete set of acenes studied here 'so!] . 
While experimental triplet energies are somewhat diffi- 
cult to compare directly with theoretical gas phase cal- 
culations, we have also included current experimental es- 
timates where available ^4^, i41, i4^, i43] [82]. 



Our DMRG calculations clearly confirm that the 
acenes maintain a singlet ground-state configuration and 
that there is a finite singlet-triplet gap for all chain- 
lengths. Going from the minimal ST0-3G to the DZ ba- 
sis and the corresponding larger double-active space, the 
singlet-triplet gap decreases by a few kcal/mol. With the 
DZ basis the hexacene DMRG gap is 17.5 kcal/mol. The 
remaining error in the DMRG calculations arises from the 
neglect of dynamical and a-n correlations, which would 
generally further decrease the gap size. However, we esti- 
mate the effect of dynamical correlation on the gap to be 
very small when using the complete (and the double com- 
plete) TT- valence space, on the order of a few kcal/mol. In 
TableUwe present additional CASSCF and CASPT2 re- 
sults (including the CASSCF and MRMP calculations of 
Kawashima et al. [25|) for the smaller acenes to estimate 
the effects of dynamical correlation. CASPT2 [4J] and 
MRMP [4^ both incorporate dynamical correlation on 
top of the CASSCF reference through second-order per- 
turbation theory. We observe in naphthalene that when 
using a complete 7r-valence space the CASSCF singlet- 
triplet gap is very close (within 1-2 kcal/mol) to the 
CASPT2 singlet-triplet gap. It is only when incomplete 
active spaces are used that the CASPT2/MRMP gap is 
significantly different from the CASSCF gap. In all cases, 
the DMRG complete and double tt- valence space gaps are 
closer to the experimental result than the MRMP gap 
in an incomplete active space. This highlights the im- 
portance of the complete 7r-valence space for 7r-electron 
excitations. 



Comparison of the UBLYP and UB3LYP gaps with 
the experimental data suggests that the DFT results are 
an underestimate. This is particularly true for UBLYP 
which substantially underestimates the gap. Surpris- 
ingly, the DFT gaps appear to increase between 10-acene 
and 12-acene. 



III. THE SINGLET-TRIPLET GAP 

In Fig. [5] we present the computed DMRG singlet- 
triplet energy gaps as a function of the acene length. 
The calculations on dodecacene correspond to a (50, 50) 
CASCI and are only made possible through the DMRG 
algorithm. Included for comparison are the UB3LYP/6- 



Using our DMRG data we can extrapolate to the in- 
finite polyacene limit. Empirically, we find that the 
singlet-triplet gap is well fitted by an exponential form 
a + be^'^, giving a gap for the infinite chain of 8.69 ± 
0.95 (ST0-3G) and 3.33 ± 0.39 (DZ) kcal/mol respec- 
tively, somewhat lower than the previous estimate of 12.2 
kcal/mol obtained by Raghu et al. 30] using the semi- 
empirical Pariser-Parr-Pople (PPP) Hamiltonian [46l.[47j. 
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TABLE I: Singlet-triplet gap (iJtripict — i^singict) energies in kcal/mol for the acene series. 
[n]acene DMRG/STQ-3G DMRG/DZ UB3LYP/6-31G(d) UBLYP/6-31G(d) Experiment 
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FIG. 2: Singlet-triplet energy gap as a function of the acene 
length. 



IV. POLYRADICAL CHARACTER OF THE 
GROUND-STATE 

Having established that the acene ground-states are 
singlets, are they then singlet diradicals as argued by 
Bendikov et al. [1^? A simple way to establish whether 
there are unpaired electrons in a correlated wavefunction 
is to examine the occupation numbers of the (spinless) 
natural orbitals - in a closed shell configuration, these 
are always 2 (doubly occupied) or (unoccupied), while 
values close to 1 indicate single occupancy and unpaired 
electrons ^Sj. In Fig. [3] we plot the occupancies of 
the natural orbitals for the acene series. We have des- 
ignated the two orbitals with occupancies closest to 1 
the HONO ("highest occupied natural orbital" with oc- 
cupancy greater than 1) and LUNO ("lowest unoccupied 
natural orbital" with occupancy less than 1) respectively. 
These natural orbitals together with usual HOMO and 
LUMO are shown in Fig. 01 

As can be seen, as we proceed to longer acenes the 



TABLE II: Effect of active-space size and dynamical correla- 
tion on the singlet-triplet gap in smaller acenes. Complete = 
complete yr-valence space, double = double 7r-valence space, 
partial = incomplete active space: 2-acene (8,8), 3- and 4- 
acenes (12,12). DZP = Dunning DZ basis with polarization 
functions [37l. [3^ except for results of Kawashima et al. [25| . 
All energies in kcal/mol. 
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45.9 


34.7 


double/DZ 
DMRG 


61.0 


44.0 
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61.0 


43.1 


29.3 



"CASSCF/MRMP calculations of Kawashima et al. [23; vertical 
singlet-triplet gap in a cc-pVDZ basis without polarization func- 
tions on H. 



occupancy of the HONO and LUNO indeed approach 1, 
which is consistent with Bendikov et aZ.'s prediction of di- 
radical character. The DZ basis, while yielding less rad- 
ical character (e.g. the occupancy of the HONO in pen- 
tacene is 1.66 and 1.73 using the ST0-3G and DZ bases 
respectively) shows the same general behaviour (The de- 
creased radical character in the DZ basis is consistent 
with the general observation that radical character is re- 
duced by dynamic correlation). However, what is surpris- 
ing is that if we follow the trend for the next nearest single 
occupancy orbitals (the HONO-1 and the LUNO-l-1), the 
rate at which they approach single occupancy is compa- 
rable to that of the HONO and LUNO. This suggests that 
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natural orbital index 



FIG. 3: Natural orbital occupation numbers for the acene 
series in the STO-3G basis. The lines are guides for the eye 
to show the evolution of the occupation numbers for the near- 
singly occupied orbitals as a function of chain length. 

if we were to proceed to acenes longer than the 12-acene, 
we would eventually find not a diradical ground-state, but 
a polyradical ground-state. 

Several different measures of the number of "effectively 
unpaired" electrons in a molecule have previously been 
proposed. While such integrated measures must contain 
less information than the underlying distribution of natu- 
ral orbital occupations examined above, we include them 
here for completeness. We have investigated two mea- 
sures, due to Takatsuka liol . Isol . Isil Is^ . iSS] and Head- 
Gordon respectively, [13, l56l |. defined via 

-^uiifSbod^'' = X^i — nf (1) 
A^^npat" =E.niin(n„2-n,) (2) 

Here Ui is the occupation number of the ith natural or- 
bital, which ranges from to 2. The contribution from 
each orbital is a maximum when Ui — 1, whence each 
orbital contributes 1 electron to the effective number of 
unpaired electrons. These measures are plotted for the 
acenes in Fig. [5j 

Some care must be taken when interpreting Fig. [S] 
Certainly, the values must not be taken literally; 12-acene 
does not contain 14 unpaired electrons! Both measures 
are extensive, meaning that they increase with the size 
of the molecule. (This also means that a large enough 
assembly of nearly closed-shell molecules would appear 
to have a substantial number of unpaired electrons us- 
ing these measures. In such a case, however, the HONO 
and LUNO occupation numbers would not change and 
would stay near 2 and as the number of molecules is 
increased, unlike what we see in the acenes). However, 
extensive scaling should only be observed for system sizes 
larger than the typical size of an unpaired electron. Ex- 
amination of the gradients of the plots in Fig. [5] shows 



an onset of extensive scaling around hexacene, which is 
consistent with the observation of Bendikov et al. that 
the first symmetry breaking of the density functional cal- 
culations occurs also at this chain-length. This further 
suggests that we can roughly associate one unpaired spin 
with every five to six rings. 

V. VISUALISING ELECTRON CORRELATIONS 
A. Theoretical background 

Let us then regard the ground-state of the longer 
acenes as being a singlet polyradical. How are we to 
understand its electronic structure? We can visualise 
the simultaneous behaviour of the multiple electrons in- 
volved in the polyradical behaviour through their correla- 
tion functions. Correlation functions have long been used 
to understand bond-alternation and electron correlation 
in conjugated systems [1^, [l^l- We have evaluated three 
correlation functions, the particle-particle, spin-spin, and 
singlet diradical correlation functions which we now de- 
scribe. In this section we shall be concerned with the 
correlations of the electrons in real-space. Thus in the 
following, indices i and j always refer to the (orthogo- 
nalised) Pz atomic orbitals on atoms i and atom j respec- 
tively, (nf) and (S^) refer to the average a occupancy 
and z-component of the spin in these orbitals. 

Particle-particle : 

Cparticic(*, j) = 4 - )) (3) 

This measures the correlation between the a population 
of orbital i and /3 population of orbital j. In a single- 
determinant wavefunction (such as the Kohn-Sham wave- 
function) there are no a/3 correlations and this quantity 
is identically zero. 
Spin-spin: 

Cspi„(*,j) = 4((5f5|)-(5f)(5|)) (4) 

This measures the correlation between the spin in orbital 
i and the spin in orbital j. Note that in wavefunctions 
that preserve the correct singlet-spin symmetry as used 
in this work, (5*1^) = (Sj) = 0. Because there are aa and 
(3P correlations from the Pauli principle even in single 
determinant wavefunctions, this quantity does not fully 
vanish in non-interacting systems. 
Singlet diradical: 

Cdiradical(*, j) = 2 [{dy^) (dT){d^} + {dU^) - {df){d^)) 

df^nf{l-n^) (5) 

The single occupancy operator df measures the proba- 
bility that an orbital i is occupied with a spin without 
any simultaneous (3 occupancy. This and the joint dirad- 
ical probability density (dfdj) were introduced by Dutoi 
et al. [11]. The function Cdiiadicai above is obtained by 
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FIG. 4: Plots of the molecular orbitals (left; LUMO+1, LUMO, HOMO, HOMO-1) and natural orbitals (right: LUNO+1, 
LUNO, HONO, HONO-1) for pentacene. 
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FIG. 5: Measures of the number of unpaired electrons in the 
acene series. See comment in text. 



removing the independent probabilities of single occupa- 
tion (e.g. {df){dj)) from the probability density of Dutoi 
et ai, to give the correlation between single occupancies 
of orbitals i and orbital j with opposite spin. Again, 
because of Pauli type correlations in single determinant 
wavef unctions, this quantity does not fully vanish in non- 
interacting systems. 

To make the meanings of these correlation functions 
explicit, we can examine the following limiting cases for 
2-electron wavefunctions (Here (f)i and (j)2 are disjoint 
orthogonal atomic orbitals). 

1. Singlet diradical 5* = ^((/>i(/>2 + 02 0i)(a/3 — /3a). 
In this case, the above correlation functions as- 
sume their extremum values of 1 or —1. Thus we 

would find Cparticlo(l, 1) = Cparticlo(2, 2) = -1, 
C'spin(lj 1) — C'spin(2, 2) = 1, CdiradicaKlj 1) = 
Cdiradical(2, 2) — —1, and Cparticlc(l, 2) = 1, 



CspinClj 2) 



-1, and C, 



diradical 



(1,2) 



1. 



2. Closed-shell singlet * = 5^((/'i+02)(0i+02)(a/3- 
Pa). In this case all correlation functions are iden- 
tically 0, reflecting the absence of a(3 correlation. 
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3. Triplet rris = diradical ^' = ^(0i02 — 4>24'i){o-P + 
f3a) . Here all correlation function values are identi- 
cal to those for the singlet diradical wavefunction. 
(One observes that the = triplet wavefunc- 
tion differs from that of the singlet diradical only 
in the phase relationship between orbital products 
such as 4>i4>2 ^^"i 0102 ' t^u^ to distinguish the two 
one should examine averages such as {S^SJ). Note 
that ms — 1 triplet states were used in this study) . 

B. Correlation functions 

Figs. [6] and [7] show plots of the different correla- 
tion functions evaluated for the singlet ground states of 
napthalene, pentacene and dodecacene. Since the corre- 
lation functions are functions of two positions, we have 
plotted them as a function of the second position with 
the first (reference) position fixed (indicated by the boxed 
value in the figures). 

In Fig. [5] all plots have the reference position fixed at 
the centre of the lower acene strand. Examining Cparticio 
we see that it is large and negative at the reference posi- 
tion. Thus given an a electron in this orbital, there is a 
significantly decreased chance of finding a simultaneous 
(3 occupation of the orbital, or more simply, double occu- 
pancy of the atomic orbital is disfavoured. Moving one 
atom away, Cparticic is large and positive reflecting an in- 
creased chance of finding the orbital to be occupied with 
opposite spin to that at the reference position. This an- 
tiferromagnetic correlation continues further away from 
the reference position in a pattern of positive and nega- 
tive values of Cparticic, though the rapidly decreasing am- 
plitudes indicate that the correlations are short-ranged. 

Examining the spin-spin Cgpin and singlet diradical 
Cdiradicai Correlation functions yields a similarly consis- 
tent picture. Cgpin is large and positive at the reference 
position while Cdiradicai is large and negative, which both 
indicate that the orbital has a strong tendency towards 
single occupation. The neighbouring atoms further show 
strong single-occupancy, antiferromagnetic correlation 
with large negative (positive) values of Cgpin (Cdiradicai), 
and this correlation decreases rapidly further away from 
the reference position. While Cgpin and Cdiradicai do not 
identically vanish for a single determinant uncorrelated 
wavefunction, their corresponding plots for pentacene in 
Fig. [6] show that aside from a small reduced propensity 
for double occupancy at the reference position which re- 
sults from electron delocalisation, there are no significant 
antiferromagnetic correlations along the chains. (Recall 
that Cparticic is identically zero in the single determinan- 
tal wavefunction). 

Comparing the correlation functions of napthalene and 
dodecacene shown in Fig. [SI for which the reference atom 
is in both cases at an "inner" position on the strand, we 
see that there is a (slight) increase in the antiferromag- 
netic correlations as the length of the acene increases. In 
napthalene the correlation between the reference atom 



and the atom on the neighbouring strand is stronger than 
the correlation to its neighbours on the same strand, a 
situation which is reversed in the longer acenes. This is 
consistent with the increasing difference between the lad- 
der and rung C-C bond-lengths, which leads to the view 
of the longer acenes as a pair of coupled polyacetylene 
strands 

Fig. [7] shows the spin-spin correlation plots where we 
move the reference position around the ring. As expected 
the antiferromagnetic correlations persist as the reference 
position is moved. Bond alternation is stronger near the 
edges of the pentacene ring and this leads to asymmet- 
rical correlations between the reference position and its 
neighbours; stronger correlations are observed across the 
shorter bonds. 



VI. THE NATURE OF BONDING IN THE 
ACENE POLYRADICAL STATE 

The correlation functions evaluated above present a 
dynamic picture of the electronic motion in the acenes. 
Tracking a single electron as it makes its way around the 
ring, a second electron is pulled along, antiferromagneti- 
cally coupled to the first and distributed over the nearest 
neighbour atoms. 

Short-range antiferromagnetic correlations naturally 
bring to mind resonating valence bonds [59, 60,, Jil, 63, 
|63| . |64 | . Recall that we can expand any wavefunction in 
terms of resonance structures, which may be classified 
as covalent, singly ionic, doubly ionic and so on (see Fig. 
|8]). (In this language, the resonance structures are viewed 
only as a many-body basis for expanding the wavefunc- 
tion; the molecular geometry is fixed across the different 
structures). From our correlation functions we see that 
the acene ground state is dominated by covalent reso- 
nance structures (no double occupancy of the Pz orbitals) 
with short-range spin-couplings (i.e. short-ranged anti- 
ferromagnetic correlations) . 

We should note that the use of the word resonance 
here is different from the colloquial usage where reso- 
nance structures are a simple metaphor for delocalisation. 
In terms of the resonance structures as a physical basis, 
electron delocalisation requires superposition of covalent 
and ionic structures with roughly equal weights. In the 
limit where the wavefunction is comprised only of cova- 
lent structures, we instead have an extreme localisation: 
the electrons are fixed and unpaired on each of the atomic 
sites with only a spin degree of freedom, which fluctuates 
between different kinds of spin-couplings. Wavefunctions 
which are predominantly covalent in nature can therefore 
be viewed as polyradicals as every electron is unpaired 
in a localised, isoenergetic atomic orbital. The covalent 
nature of the acene ground-state revealed by the correla- 
tion functions argues for this polyradical interpretation, 
which is consistent with the picture given earlier by the 
natural orbital occupations. 

Valence bond language in conjugated 7r-systems has 
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(a) 



C, 



C, 




' -n.oni -0-011 -0.020 -0.033 -0.050 -o.ore -0.076 -o.oso -0.033 -0.020 -0.011 -0.004 

•particle (l.ool 0.007 0.010 0.016 0.026 0.0-16 0.083 0.046 0.026 (1.016 0.010 0.007 0.004 




1 -0.007 -0.010 -0.015 -0.027 -0.056 l-0.32()l ^ -0.056 -0.027 -0.015 -0.010 -0.007 
0.004 0.010 0.020 0.035 0.066 0.191 0.191 0.066 0.035 0.020 0.010 0.001 



' 0.004 0.011 0.020 0.033 0.050 0.077 0.077 0.050 0.033 0.020 0.011 0.004 

•spin -0.004 -0.007 -o.oio -o.oie -0.027 -0.052 -0.I60 -0.052 -0.027 -11.016 -0.010 -0.007 -o.ooJ 




-0.004 -0.010 -0.020 -0.035 -0.068 -0.314 -0.314 -0.068 -0.035 -0.020 -0.040 -0.004 



' -0.004 -0.011 -0.020 -0.033 -0.050 -0.076 -0.076 -0.050 -0.033 -0.020 -0.011 -0.004 

•diradical o.ooi 0.007 o.olo 0.016 0.027 0.053 o.i82 0.053 0.027 0.016 0.010 0.007 0.004 




1 -0.007 -0.010 -0.015 -0.027 -0.057 ^ l-0.43al ^ -0.057 -0.027 -0.015 -O.OIO -0.007 
0.004 0.010 0.020 0.035 0.069 0.370 0.370 0.069 0.035 0.020 0.010 0.004 



(b) 

r" -0.000 -0.001 0.000 -0.001 -0.000 r" 0.000 0.001 0.000 0.001 0.000 

^spin -(I. (100 -0.001 -11.000 -0.000 -0.001 -0.00(1 diradical o.iioii 0.000 0.000 0.000 o.ooo 0.000 




-0.000 -0.015 -0.013 -0.013 -0.1115 -0.000 0.0011 0.008 0.005 0.005 0.009 0.000 



-0.001 -0.059 10.18.11 -0.059 -0.004 0.006 0.056 1-0.0311 0.056 0.006 

FIG. 6: (a) Particle-particle, spin-spin, and singlet diradical correlation functions evaluated for napthalene, pentacene, and 
dodecacene in the STO-3G basis, fixing the reference point at the centre of the lower acene strand (indicated by the boxed 
value). The value of the correlation function is indicated by the numbers; the size and colour of the circles give the magnitude 
and sign of the correlation function respectively, (b) Correlation functions in a non-interacting model of pentacene. 
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-0.532 -0.085 -0.030 -0.013 -0.004 -0.109 -0.387 - 0.065 -0.023 -0.008 0.024 0.072 |Tr705l 0.072 0.024 

pTfiSil ^ 0.064 0.027 0.014 0.008 0.005 0.004 BHMI „ 0.050 0.026 0.011 0.008 -0.030 -0.065 -0.312^^-0.312 -0.065 -0.030 






-0.192 -0.062 -0.030 -0.010 -0.009 -0.005 -0.002 -0.193 -0.063 -0.031 -0.016 -0.009 0.023 0.041 0.074 0.074 0.041 0.023 

0.081 0.044 0.023 0.011 0.004 0.047 0.079 0.041 0.020 0.008 -0.023 -0.064 -0.185 -0.064 -0.023 



-0.013 -0.048 -0.312 -0.249 -0.029 -0.008 -0.023 -0.065 -0.387 -0.169 -0.004 -0.013 -0.030 -0.085 -0.533 

0.014 0.026 0.058 10.6.531 _ 0.059 0.027 0.008 0.014 0.026 0.059 ^ BI659I 0.004 0.005 0.008 0.014 0.027 0.004 10.6961 






-0.010 -0.031 -0.062 -0.193 -0.063 -0.030 -0.009 -0.010 -0.031 -0.063 -0.193 -0.062 -0.005 -0.009 -0.016 -0.030 -0.062 -0.192 
0.013 0.035 0.074 0.064 0.024 0.008 0.020 0.041 0.079 0.047 0.004 0.011 0.023 0.044 0.081 

FIG. 7: Spin-spin correlation functions for pentacene as we move the reference point of the correlation function (indicated by 
the boxed value) around the ring. 




FIG. 8: Covalent (top) and ionic (bottom) resonance struc- 
tures for a conjugated system. Note that electron delocalisa- 
tion requires a combination of both covalent and ionic reso- 
nance structures with roughly equal weights. 



long been appreciated in the context of their low-lying 
states, which are conventionally classified as covalent or 
ionic depending on the main resonance structures. Typ- 
ically, covalent states appear more naturally described 
in the valence bond language. For example, in the 
polyenes, which may be thought of as making the two 
legs of the acene ladder, the lowest excitation is a covalent 
2Ag state which appears to have large double excitation 
character from a molecular orbital viewpoint. However, 
this low-lying double excitation is easily understood in 
the valence bond language as arising from the singlet- 
recoupling of two sing let— >triplet excitations on adjacent 
double bonds [H [gi, HE IH, [H, |l3 . Valence bond de- 
scriptions and analyses have also been examined in the 
context of radical electronic structure |18„ Jl]. In tradi- 
tional CAS calculations, to extract a valence bond picture 
one usually re-parametrises the wavefunction through a 
valence bond ex pan sion (sometimes known as CASVB 
[Zi [Zl, H III, llil)- Such studies also find that the 
benzene ground-state should be viewed as a covalent 
state with antiferromagnetic spin-couplings, in accor- 
dance with what we have found for the acenes. The expo- 
nential size of the valence bond basis limits the CASVB 



analysis to small molecules, but as we have demonstrated, 
correlation functions provide an alternative mechanism 
to infer the resonance nature of a state. 

In a simple view of bonding, such as that afforded by 
the Hubbard [t^, [ill or Pariser-Parr-Pople models there 
are two scales of energy, the resonance or hopping en- 
ergy t associated with delocalisation and the Coulomb 
repulsion energy U associated with double occupancy of 
an atomic orbital. When U/t 1, we may be said to 
be in the strongly interacting regime. Under such cir- 
cumstances, the molecular orbital picture begins to break 
down and instead the appropriate qualitative wavefunc- 
tion is the superposition of covalent resonance structures 
as described above. A standard choice of parameters 
for conjugated polymers in the Ohno parametrisation of 
the PPP Hamiltonian is U = 11.26eV and t = 2.4eV 
(29I [tqI . [soj . placing systems such as the acenes in the 
moderately strongly interacting spectrum of Hamiltoni- 
ans and therein lies an understanding of the polyrad- 
ical character and covalent ground-state that we have 
observed. 



VII. CONCLUSIONS 

In summary, motivated by predictions of unusual 
ground-states in the longer acene molecules, we inves- 
tigated acene electronic structure with high-level wave- 
function theory. Using a new ab-initio Density Ma- 
trix Renormalization Group algorithm we could carry 
out Complete Active Space calculations on the acenes 
from napthalene to dodecacene that correlated the full 
TT- valence space. We find that the ground-state remains a 
singlet as the chain-length increases, with a finite singlet- 
triplet gap in the infinite chain limit. Detailed exami- 
nation of the wavefunctions, natural orbitals, and effec- 
tive number of unpaired electrons further reveals that 
the longer acenes exhibit singlet polyradical character in 
their ground-state. Through a series of correlation func- 
tions we observe that electrons are antiferromagnetically 
coupled in pairs on neighbouring atoms as they move 
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around the acene chains. These results are consistent 
with a view of the longer acenes as moderately strongly 
interacting electronic systems, for which the appropriate 
reference description is a polyradical wavefunction aris- 
ing from a resonance of predominantly covalent valence 
bond structures. We note that such a viewpoint is es- 
sential to understand the excitations of these systems. 
Finally, our study illustrates that even simple systems 
such as the acencs can provide unusual surprises in their 
electronic structure. 
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